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Abstract 

A closed analytic expression is given for the spectrum of low energy photons in the 
annihilation of orthopositronium, which expression sums all the effects of the Coulomb 
interaction between the electron and the positron. The applicability of the formula 
is limited only by the condition w 7 <C m e . In the region w 7 3> m e a 2 the Coulomb 
interaction term gives the leading at low energy one-loop correction, proportional to 



aym e /w 7 , to the decay spectrum. The constant in uj-y radiative term in the one-loop 
correction to the spectrum is also presented here in an analytic form. 



1 Introduction 



The interplay of binding and radiative effects in the properties of positronium makes it an 
interesting case study in pure QED. In particular, the annihilation of the lowest 3 Si state 
of orthopositronium (o-Ps) into three photons has been studied both experimentally and 
theoretically since the first calculation of this process by Ore and Powell 1J. Most recently 
Manohar and Ruiz-Femema[2 have considered the effects of the Coulomb interaction on the 
photon spectrum in the o-Ps decay in the region of small photon energies, where these effects 
are essential. They have found in an integral form the expression for this spectrum, which 
takes into account to all orders the Coulomb interaction between the electron and the positron 
in the nonrelativistic approximation, and they also argued that their result is applicable in 
the region of the photon energy, uo, satisfying the condition u <C ma, with m being the mass 
of the electron. The present paper contains essentially two further developments over the 
results of Ref. [2|: a closed analytic formula in terms of just one standard hypergeometric 
function is given for the spectrum of low energy photons, and it is shown that the condition 
for its applicability is in fact uo <C m, i.e. the formula is valid in a significantly broader 
range of energies than stated in Ref. (2]. Furthermore, in the region u 3> ma 2 , the expansion 
parameter for the Coulomb effects is a^Jm/u, and the first term of this expansion describes 
the low-energy asymptotic behavior of the one-loop QED corrections to the spectrum. Due 
to the rather slow rise of the Coulomb term at low energy it is interesting to also evaluate 
the constant in uo term at small uo in the one-loop correction. As will be shown here this 
term can in fact be found from the known in the literature radiative QED corrections to the 
total rates of decay of the C even 1 Sq and 3 Po,2 states of positronium into two photons. A 
comparison of the resulting expression for the first two terms of the expansion in u/m of 
the 0(a) correction to the spectrum in the orthopositronium decay with the only available 
results of a numerical calculation jHj of this correction is also discussed here. 

2 Multipole expansion at w<m 

The spectrum of low-energy photons at w < m lends itself to a nonrelativistic treatment 1 . 

Indeed, after emission of such soft photon from the initial 3 Si state of o-Ps at energy E = 

*A similar treatment in QCD was used0] in discussion of soft gluons in a three-gluon annihilation of 
heavy quarkonium. 
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—ma 2 / '4 (relative to the threshold), the (virtual) e + e~ pair remains nonrelativistic and is in 
a C-even state and has negative energy E = E Q — u, which state then annihilates into two 
hard photons. This picture can be represented jH |2] by the diagram (of the nonrelativistic 
perturbation theory) shown in Fig.l. The annihilation into two hard gluons takes place 
at distances 0(m _1 ), which are considered as infinitesimal in terms of the nonrelativistic 
relative motion of the electron and the positron. The Green's function at energy E, describing 
the propagation of the pair between the emission of the soft photon and annihilation, thus 
contains the exponential factor exp(— nr) with k defined as —n 2 /m = E, so that the discussed 
process is determined by distances of order k~ x and the typical velocities of the electron and 
positron are given by v ~ K,/m. It is well known that in such situation the expansion 
parameter for the Coulomb interaction is am/ n. Thus in the region where k is comparable 
with ma one has to use the exact Green's function in the Coulomb potential. In the region 
of larger k: k 3> ma, the nonrelativistic treatment is still applicable as long as n 2 <C rn 2 , 
and the Coulomb effects can be calculated either by a perturbative expansion of the Green's 
function, or by an expansion of the exact formula, if such formula is available. 



I 



o-Ps 




Figure 1: The diagram for the description of the photon spectrum in the region u <C m 
in the decay of o-Ps. The open circle stands for the interaction of a nonrelativistic e + e~ 
pair with a soft photon (dashed line), and the filled circle shows the annihilation into two 
hard photons (wavy lines) at distances 0(m _1 ). G c stands for the Green's function in the 
Coulomb field. 

The Hamiltonian for the emission of the soft photon in the diagram of Fig.l can be 
expanded in multipoles. The dominance of the lowest multipoles is guaranteed by a general 
consideration jHj for as long as the system is nonrelativistic, i.e. uo <C m. However, since 
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at this point the present paper differs from Ref. [2], it is appropriate to provide here a more 
detailed discussion. The C-even state in Fig.l annihilating into two hard gluons can be either 
a spin-triplet state with an odd orbital momentum L — 2n + 1, or a spin-singlet state with 
an even orbital momentum L = 2n, where in both cases n is a non-negative integer. For the 
first set of states the minimal multipole contributing to the transition 3 S ^ 3 (2n + 1) + 7 
is of the electric type: E(2n + 1), and the amplitude of the transition is proportional to 
(ujr)( 2n+1 \ while for the latter states the lowest multipole is of the magnetic type: M(2n), 
and the amplitude of the transition is proportional to (u/m) (ujr)( 2n \ The amplitude of the 
annihilation of a state with orbital momentum L at distances 0(m _1 ) contains L derivatives 
of the wave function at the origin, i.e. it is proportional to (n/m) 1 . Multiplying the indicated 
factors for both sets of the intermediate states, and taking into account that the exp(— nr) 
behavior of the Green's function constrains the product nr at order one, one readily finds in 
both cases that the contribution of the corresponding intermediate state in the amplitude 
described by Fig.l contains the factor (u / m) ( - 2n+1 \ Thus in the nonrelativistic region of 
l) <C m it is sufficient to consider only the intermediate states with the lowest n, i.e. n = 0. 
Clearly, these states are the 1 S and 3 P, and they are reached from the initial 3 S state by 
respectively Ml and El radiative transitions. It is important to emphasize that both these 
intermediate states provide contribution of the same order in the nonrelativistic limit [HE], 
and also to notice that the spatial extent, ~ (ma) -1 , of the initial state of positronium does 
not enter as a parameter in the discussed multipole expansion 2 



3 The Coulomb interaction effect in the spectrum 

Using the standard expressions for the Hamiltonian of the Ml and El interaction, and 
also using the well known amplitudes of the two-photon annihilation of the 1 S jH] and 3 P [7] 
states, one can write the expression for the differential rate of the o-Ps annihilation in the 
limit x = uj jm <C 1 in the form[2] 



dr ma 6 
dx 9n 



1 I 2 J 


7 ! |2l 


I \ 
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2 The distances in the radiative transition amplitude are constrained by the falloff of the Green's function, 
rather than by the falloff of the wave function of the initial state. In particular at k> ma the initial state 
wave function enters only through its value at the origin. 
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where a m and a e are respectively the magnetic and electric dipole amplitudes, which can 
be written in terms of the wave function and of the discussed Coulomb Green's function 
G c (x, y; —K 2 /m) as follows, 

a m = J G c (0, y; -K 2 /m) ^ (y) d 3 y , (2) 

<My)dV (3) 



3Vo(0) 



fji 



d 

— G c (x,y; -K 2 /m) 

OXi 



x=0 



Here ipo(y) stands for the wave function of the initial state of the orthopositronium, and 
the factor i/jq(0) appears in the denominator in both these formulas due to that its value 
is already included in the normalization in eq.(fTj). In other words, the amplitudes a m and 
a e are normalized in such a way that they both are equal to one if the Coulomb Green's 
function G c is replaced by the free motion one Gf, in which limit the lowest-order formula[l 
for the spectrum is reproduced. 

The magnetic amplitude a m is however trivial and is equal to one also if the Coulomb 
interaction is taken into account. This clearly is a consequence of that the spatial wave 
function of the ground 3 S state is orthogonal to those of all the 1 S states except for the 

ground one, where the overlap integral is equal to one. The relation CLjj-y X IS valid up to 

relativistic terms, including the 3 S — 1 S hyperfme splitting, which effects are beyond the 
intended accuracy 3 . 

The integral in eq.(j2J can be calculated using the partial wave expansion of the Green's 
function 

G c (x, y; E) = ]T(2£ + 1) G e (x, y; E)P £ ( ^ ] , (4) 

with Pg(z) being the Legendre polynomials, and the following representation [8J of the partial 
wave Green's functions Ge in the Coulomb problem: 

G e x,v; = (2kxY (2KyYe K(x+y > V ; n , v J" - , (5) 

where v = ma /(2k), and L^\z) are the Laguerre polynomials defined as 



3 The hyperfine splitting, becoming essential at very small u>: uj ~ ma 4 , was considered in Ref.[2] in order 
to establish the low-energy behavior mandated by the Low theorem. 
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(Usage of this representation can also be found in similar problems in Refs. 0] and 0.) 

The expression in eq.(J3} contains only the P wave partial Green's function G\ and thus 
can be written as 

47T uj r°° ( k 2 \ 



■00 (0) jo \ ' ' m 



24 ^ n + 2 - z/ Jo "V 2 



0, - 



n + 4)!/ l-z/\ n l + z/(n + 3)!/ 1 - i/ 



24n! V 1 + vJ 2 6n! V 1 + v 



(7) 



Here a use is made of the explicit form of the ground state wave function: ipo(y) I Vo(0) = 
exp(— may/2) as well as of the relation mu/n 2 = 1 — v 2 . In the last transition the formula 
(|)J) is used to perform integration by parts. The sum in the latter expression in eq.([7]) is of 
the Gauss hypergeometric type and can be done explicitly in terms of the hypergeometric 
function 2 iq, so that the final result can be written as 

, N (l-v)(3 + Bv) 8z/ 2 (l-z/) / \ — v 



3(1 + v) 2 3 (2 - z/) (1 + z/) 3 z 1 V '' ' l + i/. 
The explicit relation between the photon energy uj and the Coulomb parameter v reads as 

ma 2 1 — v 2 1 1 — v 2 

The amplitude a e given by the formula © is fully equivalent to its integral representation 
described in Ref. j2]. 



4 One-loop QED correction to the photon spectrum at 

The result in eq.(|Hl) can be expanded in v. The linear term in this expansion describes the 
correction of the first order in a: 

a e = l-^ + 0(u 2 ) = l-^M+0(a 2 ) . (10) 
3 3 V uj 

When used in eq.Q for the differential decay rate the linear term in this expansion deter- 
mines the asymptotic behavior of the one-loop QED correction to the spectrum at w <C m. 
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However, the slow \j\fx rise of this correction toward small x makes potentially important 
the higher terms of the expansion in x of the 0(a) one-loop correction. Clearly, the rela- 
tivistic expansion for the discussed terms due to the Coulomb interaction goes in powers of 
k 2 /m 2 « ijjjm = x, so that the next term of this type is proportional to oty/x. The same 
behavior is true for the correction terms arising from the full relativistic scattering kernel 
for the initial state and from the process (o-Ps) — > 7* — > 37. Although the latter two mech- 
anisms together account for almost 85% of the 0(a) correction to the total ratejlOj. in the 
spectrum at low x their relative contribution starts only as a^fx. Such terms are beyond 
the scope of this paper. 

On the other hand, the constant in x term, i.e. of order ax° at small x, can be quite readily 
deduced from the known in the literature results for the one-loop corrections to the total 
rates of the two photon annihilation of the C-even S and P states of the orthopositronium. 
These corrections are of genuinely radiative nature and arise from distances of order 1/m, 
as opposed to the so far discussed Coulomb effects determined by the electron-positron 
interaction at distances of order 1/k. The sources of such correction terms can be easily 
identified from the graph of Fig.l. Indeed, the elements of the considered process determined 
by the distances of order 1/m are the amplitude of the annihilation into two hard photons and 
the electron magnetic moment entering the Hamiltonian of the Ml interaction. Accordingly, 
in order to include the discussed radiative corrections, the equation (JJJ should be rewritten 
as 



dr ma 6 / , 2 fg e \ 2 r(% ^ 2 7 ) 

-x { \a m \ — — — -—+ a e 



dx 9tt 1' 1 \2/ r ( 1 SW2 7 ) 



r( 3 p -> 27) + 4 r( 3 p 2 - 27) 



r ( 3 P -> 27) 3 r ( 3 P 2 - 27) 



(ii) 

where g e is the electron gyromagnetic ratio, and the ratia of the decay rates for the indicated 
C-even states are to their values (ro) in the lowest order 4 . The 0(a) terms in these ratia 
are known in the literature: 



r(% -> 27) 

r ( l S -> 2 7 ) 
for the parapositronium decayjTT], and 




(12) 



r( 3 P ^2 7 ) a /tt 2 7\ T( 3 P 2 -> 2 7 ) Aa 

r ( 3 Po^27) % 4 3 r r ( 3 P 2 ^2 7 )^ tt 1 ] 



4 Thcre is of course no contribution from the 3 Pi state due to the Landau- Yang theorem. 
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for the 3 -Po,2 states, which formulas are an adaptation to QED of the quarkonium results for 
QCD corrections in Ref . 1 1 2| . 

Using the expressions (|T2*|) and (fT3~|) and also the famous Schwinger's result g e /2 = 1 + 
a/ (27r)+0(a 2 ) in eq.fjllj). the formula for the spectrum can be written including the radiative 
correction: 



dr 

dx 



ma 



9tt 
10 ma 6 
27tt 



-x 



1 + 



a 



7T 



-X 



71 



14tt 



7T Z 
T 

4 3tt 2 



-4 



+ ke 



7 a 

3 71 



71 

T 



23 

y 



15 y/x 20 



5 + o<a> 



(14) 



where the last expression also includes the first Coulomb correction to a e from eq.(|5|). 

A general calculation of the one-loop QED correction to the spectrum in o-Ps decay has 
been done numerically by AdkinsjHj, and it is instructive to compare the two results. For 
the purpose of such comparison, following the conventions of Ref.jS], we write the formula 
for the differential probability in the form 



jL_dT 

r dx 



a 



7o(x) +-7i(x) . 

71 



(15) 



Here T = 2 (tt 2 — 9) ma 6 / (9tt) is the total decay rate in the lowest order pQ, 7 is the normal- 
ized differential decay rate in the same lowest order, for which we use here its nonrelativistic 
limit at small x: 7o(x) = (5/3)x/(7r 2 — 9), and ji(x) is the first-order QED correction to the 
spectrum. From the equation (|T4*j) one readily finds that under the described conventions 
7i(x) is given by 



7i0) 



9 (tt 2 - 9) 





3tt 2 \ 1 


U7Cy/x + 15 ^ - 


X 




20 J 



(16) 



In Ref. |S] the physical region of x from x = to x = 1 is divided into 20 equal bins of 0.05 
each, and the integral of 71 (x) over each of the bins is tabulated. In particular in the first 
two bins the integrals are found as respectively -0.0502(11) and -0.0994(14). The formula in 
eq. (|16[) gives for the same integrals numerical values of -0.0467 and -0.0911. The difference 
from the result of the full calculation j^j, albeit numerical, is in a reasonable agreement 
with the accuracy expected from using only the first two terms of the expansion in y/x at 
x w 0.05 - 0.1. 
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